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Abstract 

\Q An introduction and summary is given of the main results achieved by working group 1: 

Structure Functions in Deep Inelastic Scattering at HERA. The prospects were discussed of 

^— l future measurements of the structure functions F2, Fl, xGs, F®® and at HERA. 

O | The results represent a long term programme of experimentation with high luminosity, 

D different lepton beam charges, proton and deuteron beams, allowing for precision mea- 

surements. The theoretical investigations focussed on QED corrections, higher order QCD 
corrections for different observables, resummation of small x contributions, and the de- 
tailed understanding of NLO QCD evolution codes, to allow for a high precision analysis 
of the forthcoming deep inelastic data. 

O 

\o 
o\ 

^ 1 Introduction 
1 

<D The discussion of prospects of structure function physics has taken place for the third time in a 
HERA physics workshop. In 1987, the emphasis was put mainly on high Q 2 structure function 
simulations, and the determination of part on distributions [1]. In 1991, the experimental part 
concentrated on the reconstruction of the deeply inelastic event kinematics and Monte Carlo 
simulations, while the theoretical investigations began to focus on the physics of the parton 
densities at very low x [2] . 

This workshop was based on the experience gained in the first four years of running at 
HERA, and on recent important theoretical developments, in particular the extension of next- 
to-leading order (NLO) calculations to more observables, and on small x resummations. One of 
the major goals on the theoretical side was to summarize the progress towards high-precision 
calculations, and to point out which further theoretical investigations are most important in 
order to match the accuracy expected in future high-statistics structure function data from 
HERA. 
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1 Summary report of Working Group 1 of the 1996 HERA Physics Workshop 



On the experimental side, the scope of the studies has been extended to the charm con- 
tribution to the proton structure functions at low a:, and to the pion structure function. The 
knowledge of existing HERA data allowed for much more reliable simulations than in previous 
workshops. Demands on the future HERA running were derived, which are contained in the 
experimental section of this paper. A challenging scenario of running HERA for many years has 
been developed. It requires the accumulation of high luminosity for both lepton beam charges 
and also to accelerate deuterons. This will permit to reach very high precision, i.e. an error 
level of only a few % for the proton stucture function F 2 , and also accurate results for Fl, F% 
and Fi,-. On the basis of this data high precision studies of QCD in the regions of very small 
x and of very large Q 2 can be carried out. In the following, we summarize the main results 
obtained in the theoretical and experimental investigations of the working group. 

2 Theoretical Studies 

Since the 1991 HERA Physics Workshop, various theoretical developments have taken place 
in the field of deep-inelastic scattering (DIS). These include the extension of QED radiative 
corrections to more observables, the calculation of various QCD coefficient functions and scat- 
tering cross sections to 0(a 2 ), determinations of the 3-loop corrections of a series of moments 
of coefficient functions and anomalous dimensions for non-singlet as well as singlet structure 
functions, and the resummation of certain classes of logarithmic corrections at small x. Dur- 
ing the present workshop, applications of these results have been considered to future HERA 
data. Comparisons have been performed of the experimental requirements and the theoretical 
accuracy achieved at present. 

2.1 QED Radiative Corrections 

The understanding of the 0(a) QED corrections for a few choices of kinematical variables was 
already obtained at the time of the 1991 HERA workshop [3]. In a contribution to the present 
workshop [4], recent developments in the semianalytical approach by the HECTOR collabora- 
tion [5] have been summarized. These include the 0(a) corrections for a wide variety of kine- 
matical variables as well as leading-log 0(a 2 ) results. Recently also the QED corrections due 
to the resummation of the 0(a[aln 2 (l/x)] 1 log[Q 2 /m 2 ]) terms were calculated [6]. For leptonic 
variables, their contribution is of similar size as the second order corrections ~ a 2 \og(Q 2 /ml) 
and diminishes the effect of these terms. A first complete 0(a) calculation was reported also for 
neutral current polarized-lepton polarized-hadron scattering [7] accounting for photon and Z- 
boson exchange. This calculation covers the cases of both longitudinal and transverse nucleon 
polarization for the twist-2 contributions to the Born term [8]. 

2.2 Comparisons of NLO QCD Evolution Codes and the Theoretical 
Error of a s (M 2 z ) 

A major effort within the working group concerned a detailed comparison of next-to-leading 
order (NLO) QCD evolution codes. The aim of this comparison was to understand the ac- 
curacy of the different numerical solutions of the evolution equations and their conceptual 
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differences. A numerical agreement of better than ±0.05% of the parton densities has been 
achieved between five of the evolution programs. For these codes, previous deviations due to 
different theoretical assumptions on the truncation of the perturbative series at NLO are now 
completely understood. The results of this study are summarized in ref. [9]. 

Precision measurements of F 2 allow for detailed QCD tests and an improved extraction of 
a s (Q 2 ). The theoretical uncertainties of the a s measurement from scaling violations have been 
investigated for NLO analyses in the HERA range [10, 11]. The different representations of 
the evolution within NLO lead to a shift of up to Aa s (M§) = 0.003. The largest theoretical 
errors at NLO are due to the renormalization (R) and factorization (M) scale uncertainties, 
resulting in Aa,(Mf) = \ R \ M , for a Q 2 -cut of Q 2 > 50GeV 2 . The contribution 

of mass threshold uncertainties to Aa s (M§) was conservatively estimated to 0.002 in ref. [12]. 
Due to the large statistics at low Q 2 it appears to be desirable to include also the range down 
to a few GeV 2 in the QCD analyses. To exploit this region fully, however, requires to carry out 
next- to-next- to leading (NNLO) analyses, for which the 3-loop splitting functions still have to 
be calculated. 

2.3 0(a 2 s ) and O(c^) Corrections 

Since the 1991 HERA Physics workshop, various two-loop calculations have been performed 
for quantities related to structure functions. A survey on the status of these calculations, and 
recent three-loop results, was presented to this working group [13]. The NLO corrections turn 
out to be essential for the quantitative understanding of most of the observables. 

Numerical studies during this workshop were devoted to the behaviour of Fl(x, Q 2 ) [14] and 
F®® [15, 16]. An update of the parametrization of the NLO heavy flavour coefficient functions 
is also given in [15]. Moreover, a calculation of the heavy flavour structure functions in charged 
current interactions has been presented [17]. Besides of the twist-2 contributions to the structure 
function Fl, also higher-twist terms for its non-singlet part have been investigated, accounting 
for renormalon contributions [18]. Also the NLO corrections to the scattering cross section 
a (l + 9 ~ > J/^P + ^0 are available now [19]. 

First results on the behaviour of structure functions in NNLO were reported by the NIKHEF 
group [20, 21]. The first moments for both the non-singlet and singlet combinations of unpo- 
larized DIS structure functions have been calculated. By the same group also a phenomeno- 
logical analysis was presented [22], estimating the ^-dependence of the corresponding splitting 
functions by a fit allowing for a set of functions which are known to contribute. Also a phe- 
nomenological application of the non-singlet results to xF 3 was reported [23]. 

2.4 Resummations for Small x 

The measurement of DIS structure functions remains one of the major methods to investigate 
the small- a; physics at HERA. Various aspects have been considered in the working group. 
A possible approach has been discussed which relies on the k± -factorization method [24]. It 
consists of combining systematically the resummation of the small- a; logarithmic corrections, 
as given by the BFKL formalism [25] and the QCD corrections to it, with the QCD mass 
factorization theorem, dictated by the renormalization group. This approach enables one to 
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study the small- a; effects by solving improved evolution equations which include resummed 
kernels. A first numerical analysis along these lines was carried out in ref. [26]. A review of 
these equations and the current status of resummed calculations, covering also the non-singlet 
cases, can be found in ref. [27]. 

During the workshop numerical studies of structure functions at small x have been performed 
by two groups [28, 29]. One contribution to the workshop [30] dealt with the resummation of 
the small- a; contributions on the level of a double-log approximation. In the flavour non-singlet 
sector resummation effects are small, less than 1% [31]. In the singlet sector, however, they may 
give rise to large contributions [28, 29]. The question of assessing the importance of unknown 
sub-leading terms has also been addressed in ref. [28], comparing several different models. The 
outcome is that these terms seem to be able to affect the result sizeably. This indicates that 
at present the uncertainties on the theoretical predictions at small- a; are fairly large, and more 
accurate calculations (next-to-leading small- a; logarithms as well as exact 3-loop contributions) 
are necessary. The analysis in ref. [29] has emphasized the role of a combined determination of 
F 2 and Fl to pin down the behaviour of the QCD perturbative series at small x. The effects 
of small- a; resummations on the photon structure functions were also discussed [28]. 

On the leading-order level, predictions for the structure functions, covering the BFKL effects 
at small a;, can also be obtained starting from a different equation, which was introduced a few 
years ago [32] to describe the detailed structure of the gluon radiation associated to small- a; 
events. In ref. [33], the solution to this equation has been investigated numerically for the 
structure function F 2 . 

In the working group also the properties of the BFKL resummation equation itself were dis- 
cussed with emphasis on the transverse momentum cut-off [34] . Recent progress in calculating 
NLO corrections to the BFKL kernel was reported in ref. [35], see also [27] for other ongoing 
investigations. 

An alternative formulation of QCD at small x has been proposed in a series of papers, based 
on a colour dipole concept [36]. The theory has been worked out for the case of scattering 
of two quarkonia, in which non-perturbative effects are suppressed by the smallness of the 
quarkonium radius. First attempts to apply the colour dipole formulation to deep-inelastic 
scattering have been reported at this workshop [37, 38]. The suitability of this approach to 
investigate unitarity corrections and parton saturation has been emphasized, and an explicit 
parametrization of multi-pomeron exchange contributions has been presented [38] . 

3 Experimental Studies 

The future HERA measurements of structure functions using electron and positron, proton and 
deuteron beams promise to be of great interest, since the envisaged increase of luminosity will 
allow for reducing the statistical and the systematic errors, especially for the proton case, to the 
level of a few per cent in almost the full accessible kinematic range. Within this working group, 
detailed simulation studies have been performed of various structure function measurements, 
in order to estimate their expected accuracy and to analyze their physics impact. 
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3.1 The Proton Structure Function F 2 



A thorough simulation has been carried out [39] of future HERA measurements of the proton 
structure function F 2 (x,Q 2 ) for a nominal kinematic range given by y < 0.8, e < 177°, 
Q 2 > 0.5 GeV 2 , and Oh > 8°. The following assumptions were made on the future measurement 
accuracies: 0.5-1% for the scattered electron energy E' e , 0.5-1 mrad for the polar angle # e , 2% on 
the hadronic energy Eh, 1-2% for the photoproduction background uncertainty at high y, and 
2% due to trigger and detector efficiencies. The luminosity was assumed to be known within 
1%. Moreover, control of the radiative corrections at the level of 1% has been assumed. These 
conditions lead to an estimated systematic error of F 2 of about 3% in almost the full kinematic 
range of 2 • 10~ 5 < x < 0.7 and 0.5 < Q 2 < 5 • 10 4 GeV 2 . The anticipated accuracy represents 
an improvement by about a factor of two compared to the present results but extended to a 
much wider region. The current HI and ZEUS structure function analyses served as a basis for 
the simulations and were summarized in various talks presented to this working group [40]. 

In a detailed NLO QCD analysis [39], possible determinations of the strong coupling con- 
stant a s and of the gluon distribution have been considered using the HI and the ZEUS QCD 
analysis programs and fitting techniques. The error of a s strongly depends on the minimum 
Q 2 which can be included into such analysis. While perturbative QCD seems to work down 
to 1 GeV 2 at low x, the theoretical scale uncertainties in NLO become very large, see above. 
With only HERA proton data an a s error of about 0.004 can be expected for Q 2 > 3 GeV 2 . 
As described in [39] the measurement of a s requires very accurate control of the dependence 
of the systematic errors on the kinematic variables. This will permit to largely reduce their 
effect in the fit procedure like it has been practiced in the BCDMS/SLAC data analysis [41] 
which lead to an experimental a s error of 0.003. Combination of the high x fixed target F 2 
data with the low x HERA data promises to yield a precision measurement of a s with an esti- 
mated experimental error of 0.0013 for a s (M§). Simultaneously the gluon distribution can be 
measured very accurately with an estimated error of e.g. 3% for x = 10~ 4 and Q 2 = 20 GeV 2 
[39] using HERA data only. This potential measurement can only be reliably interpreted if the 
theoretical description is extended to NNLO. 

3.2 High Q 2 Structure Functions 

For Q 2 > 500 GeV 2 and with both electron and positrons beams employed, there are two 
neutral current and two charged current reactions, i.e. four cross sections to be measured. 
This permits to extract various structure functions and combinations of parton densities, as 
was already thoroughly studied in the previous HERA workshops [42, 43] and summarized in 
the discussion [44]. The main conclusions remain valid: i) the charge asymmetry in neutral 
current e ± p scattering allows to determine the interference structure function xGs = 2x ^ Q q a> q - 
(q — q) for large y; ii) the charged current cross sections directly determine the valence quark 
distributions for x > 0.3; iii) with a high precision measurement of the four cross sections, 
various combinations of parton distributions can be unfolded, e.g. the singlet distribution or 
the strange quark density. All these measurements require highest luminosities, C > 500 pb -1 , 
about equally shared between electron and positron runs. 
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3.3 Deuteron Structure Functions 



Interest has been expressed in the structure function working group [39, 45] in measuring the 
electron-deuteron DIS cross section at HERA. At low a:, despite a few % structure function 
extraction uncertainty due to shadowing effects [46], the proton-deuteron F 2 difference will 
permit to constrain the up-down quark difference. This is expected to vanish towards low x 
but has not been measured yet in the domain of x < 10~ 3 . Deuteron structure function data 
are important for a self consistent QCD analysis of HERA data as they constrain the non- 
singlet distributions in the a s and gluon determination. Even with a luminosity of C = 50pb -1 
only, interesting parton distribution combinations as s — c [43] at high Q 2 > 100 GeV 2 can be 
measured. 

3.4 Changing the Beam Energies 

With the maximum possible energies of E p ~ 1000 GeV and E e ~ 35 GeV, the cms. energy 
squared can be increased by about a factor of 1.5 as compared to the nominal energy runs. 
Hence x values lower by this factor can be accessed at a given Q 2 . It should be stressed that this 
high-energy option is the only way to explore this extended kinematical range in the foreseeable 
future. It requires only modest luminosity C ~ 5 pb -1 for the exploration of the low x region. 

Because of the relation Q 2 ~ (2E e cot # e /2) 2 , the decrease of the electron beam energy can 
be employed to reach very low values of Q 2 < 1 GeV 2 . Using a small luminosity, C ~ 1 pb -1 , 
only this will allow for a high precision measurement of the transition region to photoproduction 
even after the foreseen luminosity upgrade which may limit the detector acceptance close to 
the beam pipe. 

A run at lowest possible proton beam energy with C ~ 10 pb -1 is required in order to ensure 
maximum overlap of the HERA F 2 data with the fixed-target results. This is important for the 
precision measurement of a s . Lowering E p is also a necessity for measuring the longitudinal 
structure function Fl- 

3.5 The Longitudinal Structure Function Fl 

Various, partially complementary studies were presented to the workshop on the estimated 
measurement accuracy of the longitudinal structure function Fl [47, 48]. The joint conclusion 
of these studies is that an Fl precision measurement requires a set of about four different 
proton beam energies, E p = 250,350,450 and 820 GeV, for instance, with luminosities around 
5-10 pb -1 per energy setting. In a wide range of Q 2 , between about 4 and 100 GeV 2 , Fl 
should be measurable with an absolute error of typically 0.08 [48]. That accuracy largely 
depends on the maximum y, i.e. the minimum electron energy which can be controled in view 
of background and trigger requirements. The availability of various proton energies and the 
use of the subtraction method [49] will permit a significant measurement of the x dependence 
of Fl(x,Q 2 ) which is important to distinguish between different predictions for Fl- A run at 
highest proton energy, E p ~ 1000 GeV, would represent an important gain in sensitivity to Fl, 
and lower E e data would provide useful systematics cross checks. 
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3.6 The Charm Contribution to F 2 



At low a:, the charm contribution to the inclusive DIS cross section has been measured to be 
about one quarter [50]. Its understanding is crucial for the interpretation of F 2 and interest- 
ing for an independent measurement of the gluon distribution. With a luminosity of several 
100 pb -1 , the charm structure function will be measurable with an estimated accuracy of 10%, 
the error being mainly due to the limited knowledge of the charm fragmentation probability 
P(c — > D) into D mesons, and the detector and analyses uncertainties [51, 15]. This will permit 
a detailed investigation of the charm production mechanism as a function of Q 2 and x and a 
complementary determination of the gluon distribution. The measurement relies on the obser- 
vation of open charm production at HERA which should profit from Silicon detectors installed 
near the interaction vertex in order to enhance the charm tagging efficiency. 

3.7 Structure of the Pion 

HERA provides an interesting opportunity to study the structure of the pion [52, 53] as was 
presented to this working group [54]. Finite values of momentum transfer squared from the 
incoming proton beam to the outgoing neutron can be measured in the HI and ZEUS neutron 
calorimeters with mean values of the order of 0.1 GeV 2 . Studies were performed prior to the 
workshop [53], and a complete Monte Carlo simulation of the performance of the HI forward 
neutron calorimeter is described in [55]. A modest luminosity of 10 pb -1 would yield about 
5000 events in the range 10 < Q 2 < 15 GeV 2 . This would enable the structure of the virtual 
pion to be determined as a function of the longitudinal momentum of the exchanged pion for 
0.6 > Xir > 10~ 4 . Higher luminosities may allow the separation of the longitudinal structure of 
the pion and the determination of the structure of higher mass resonances. 
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